Critical behaviour and scaling functions of the three-dimensional 0(6) 

model 
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We numerically investigate the three-dimensional 0(6) model on 12 3 to 120 3 lattices within the 
critical region at zero magnetic field, as well as at finite magnetic field on the critical isotherm and 
for several fixed couplings in the broken and the symmetric phase. We obtain from the Binder 
cumulant at vanishing magnetic field the critical coupling J c = 1.42865(3). The universal value of 
the Binder cumulant at this point is g r (Jc) = —1.94456(10). At the critical coupling, the critical 
exponents 7 = 1.604(6), ft = 0.425(2) and v = 0.818(5) are determined from a finite-size-scaling 
analysis. Furthermore, we verify predicted effects induced by massless Goldstone modes in the 
broken phase. The results are well described by the perturbative form of the model's equation of 
state. Our 0(6)-result is compared to the corresponding Ising, 0(2) and 0(4) scaling functions. 
Finally, we study the finite-size-scaling behaviour of the magnetisation on the pseudocritical line. 
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to universal finite-size-scaling functions. We 
have indeed found that the finite-size-scaling 
functions on the pseudocritical line of 0(2) and 
0(4) are compatible with staggered QCD data. 
It has turned out problematic to check QCD 
data for further critical behaviour found in spin 
models, e.g. the goldstone effect. In QCD with 
fcrmions in the fundamental representation the 
chiral and the deconfinement phase transitions 
occur at the same temperature. This could 
change the properties of the chiral phase tran- 
sition, as additional degrees of freedom are re- 
leased. 

However in QCD with fcrmions in the adjoint 
representation (aQCD) the two phase transi- 
tions are separated, so one can study them 
individually. Since the left-handed and right- 
handed spinors are indistinguishable in the ad- 
joint representation, the chiral symmetry group 
is SU(2N f ) and not SU{N f ) L x SU(N f ) R as in 
the fundamental representation. Thus, for two 
flavours the symmetry group is SU(4), which is 
isomorph to 5*0(6), a subgroup of 0(6). The 
universality class of QCD with adjoint fcrmions 
therefore has to be that of the 0(6) spin model. 
In this paper our results for the universal prop- 
erties of this model will be shown, especially 
the scaling functions, which are needed for our 
forthcoming study of aQCD, continuing the 
work of Karsch and Lutgemeier • 
The model we investigate is the standard 0(6)- 
invariant nonlinear er-model, which is defined 
as 



I. INTRODUCTION 

The chiral phase transition of Quantumchro- 
modynamics (QCD) is of great interest for the 
understanding of the early universe and the 
physics of heavy ion collisions. For two massless 
quark flavours it is supposed to be of second or- 
der and to lie in the same universality class as 
the three-dimensional 0(4) spin model 
If these assumptions are valid, one can use the 
knowledge of the spin model to understand the 
critical behaviour of the QCD phase transition. 
It has been shown |4J, |5( that the scaling be- 
haviour in lattice simulations with two light 
quark flavours is indeed comparable to the uni- 
versal infinite volume scaling function of the 
0(4) spin model if one uses Wilson fermions, 
although the Wilson fermion action has no chi- 
ral symmetry on the lattice. For staggered 
fcrmions however the 0(4) scaling function 
does not match the QCD data. For two flavours 
on the lattice the staggered fermion action has 
a remaining U(l) x U(l) chiral symmetry. As 
this symmetry lies in the same universality class 
as the three-dimensional 0(2) spin model, the 
data has also been compared to approximated 
infinite volume 0(2) data. But the 0(2) scal- 
ing function matches even worse than the 0(4) 
function. 

Since the lattice sizes used in QCD are rather 
small, it might be better to compare the data 
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Here x and y are the nearest-neighbour sites 
on a three-dimensional hypercubic lattice, <p x 
is a 6-component unit vector at site x and H 
is the external magnetic field. The coupling 
constant J is considered as inverse tempera- 
ture, therefore J = l/T. An additional term 
+ \(4>l - l) 2 ] is often used in the Hamil- 
tonian with A tuned to minimize leading order 
corrections to scaling. It is not applied here, 
because the appropriate A value of the 0(6) 
model has not been calculated yet, and such 
a calculation is beyond the scope of this paper. 
If H = \H\ is non-zero we can decompose the 
spin vector <f> x into a longitudinal (parallel to 
the magnetic field H ) and a transverse compo- 
nent 

$x = 4>le H + & with e H = H/H . (2) 

The order parameter of the system, the mag- 
netisation M, is the expectation value of the 
lattice average 0" of the longitudinal spin com- 
ponent 

M = <^£*i> = ( 3 ) 

X 

V = L 3 is the volume of the lattice with L 
points per direction. 

At zero magnetic field (H = 0) there is no 
special direction and the lattice average of the 
spins 

<? = ^£4 (4) 

X 

will have a vanishing expectation value on all 
finite lattices, ( <f> ) = 0. As an approximate 
order parameter for M at H = one can take 

8] 

M ^ (|0|) . (5) 

Nevertheless, we can use <j> to define the suscep- 
tibilities and the Binder cumulant by 
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In the following section we describe our simu- 
lations at zero magnetic field and estimate the 
critical coupling J c from the Binder cumulant, 
the magnetisation and the susceptibilities. In 



Section [H] the critical exponents v : (3 and 7 
are determined. With simulations at H > 
111 Section [□] we investigate the behaviour of 
the model on the critical line, in the broken 
phase and in the symmetric phase. Finally, the 
resulting data is used in Section Hvl to generate 
the infinite volume scaling function of the mag- 
netisation. Using this data the infinite volume 
scaling function of the susceptibility and the 
position of the pseudocritical line are derived 
in Section A summary and our conclusions 
are given in Section IVII 



II. SIMULATIONS AT H = 

All our simulations were done on three-di- 
mensional cubic lattices with periodic bound- 
ary conditions. We used Wolff's single clus- 
ter algorithm as we did in our previous pa- 
pers (Refs. and The H = data 
were taken from lattices with linear extensions 
L = 12,16,20,24,30,36,48,60 and 72. Be- 
tween the measurements we performed 300-600 
cluster updates to reduce the integrated auto- 
correlation time Ti nt fo r the energy. 
Butera and Comi determined the critical 
point of the 0(6) spin model using a high tem- 
perature (HT) expansion as J c = 1.42895(6). 
Therefore we generally scanned the range from 
,/ = 1.3 up to J = 1.55 on smaller lattices with 
careful regard to the critical region close to the 
J c -value found in ^1 f° r an lattices. This data 
was then further analysed using the reweighting 
method. More details of the simulations near 



L 


J-range 


Nj 


N meas [1000] 




12 


1.42830 - 1.42900 


25 


100 - 200 


< 3 


16 


1.42840 - 1.42880 


18 


100 - 200 


£ 4 


20 


1.42840 - 1.42885 


19 


100 - 200 


£ 6 


24 


1.42835 - 1.42885 


19 


100 - 200 


£ 6 


30 


1.42840 - 1.42885 


18 


100 - 200 


< 6 


36 


1.42840 - 1.42880 


17 


100 


£ 5 


48 


1.42840 - 1.42880 


17 


100 


£ 6 


60 


1.42840 - 1.42880 


16 


80 


£ 8 


72 


1.42840 - 1.42880 


16 


80 


£ 9 



TABLE I: Survey of the Monte Carlo simulations at 
H = for different lattices. Here Nj is the number 
of different couplings performed in the appropriate 
J-range; Tint is the integrated autocorrelation time 
for the energy and N meas the number of measure- 
ments per coupling in units of 1000. 
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the critical point are presented in Table Q] 



A. The Critical Point T c 

Obviously any determination of critical values 
as well as the definition of the reduced temper- 
ature 



t 



T-T c 
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relies on the exact location of the critical point. 
Since there is no result from numerical studies, 
we check the aforementioned value of Butera 
and Comi first. We determine T c by studying 
the Binder cumulant g r , which is a finite-size- 
scaling function 



= Q g {tL 1/v ,L-") 



(10) 



The function Q g depends on the thermal scal- 
ing field and on possible irrelevant scaling fields. 
In this case only the leading irrelevant scaling 
field proportional to is specified, with an 
unknown u) > 0. Therefore, at the critical point 
(t = 0) g r ought to be independent of L apart 
from corrections due to these irrelevant scal- 
ing fields. Fig. H ( a ) shows our results for g r . 
On the scale of Fig. 2] (a) we observe no devia- 
tion from the scaling hypothesis. After a blow- 
up of the close vicinity of the critical point, as 
shown in Fig. Q](b), one sees that the intersec- 
tion points Ji p between the curves of different 
lattices are not coinciding perfectly at one J. 
These minor corrections to scaling have to be 
considered. By expanding the scaling function 
Q g to lowest order in both variables one gets for 
the intersection point Jj P of two lattices with 
sizes L and L' = bL 



J ip (L,b) = J c + c\s{L,b) 



with 
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To have an unbiased estimate of J c we choose 
Binder's approximation ^| 
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which can be used without knowing the val- 
ues of v and u>. In Fig. the l/Ji P val- 
ues of the intersection points from the lattices 
L = 12,16,20,24,30,36 with all other larger 
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FIG. 1: (a) The Binder cumulant g r from Eq. (JHJ 
as a function of the coupling J. The points are con- 
nected by splines to guide the eye. With increasing 
lattice size L = 12, 16, 20, 24, 30, 36, 48, 60 and 72, 
the slope of the respective curve in the critical re- 
gion increases. The vertical dashed line denotes our 
final result for J c . (b) is an enlargement of (a) near 
the critical point. The dotted lines accompanying 
the solid/dashed lines show the jackknife error cor- 
ridors. 



lattices are plotted as a function of the variable 
l/ln& of Eq. (|13|l . Linear fits should lead to 
the same value 1/J C within the errors. Fitting 
the results to a constant value and varying the 
used L- values we find 



1 



0.699960(14), (14) 

which is equivalent to 

J c = 1.42865(3). (15) 

This result agrees in the first four digits with 
the result J c = 1.42895(6) of Butera and Comi 
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FIG. 2: The coupling J ip at the intersection point 
of g r {L) and g r (bL) for various L and 6 as a function 
of 1/ In b, Eq. 1131 . The lines are linear fits for the 
lattices L = 12 (black), 16 (red), 20 (green), 24 
(pink), 30 (dark blue) and 36 (light blue) through 
the intersection points with all larger lattices. The 
errors of the endpoints are drawn on the ordinate. 

[l2|. There is a slight difference in the last two 
digits. As this difference is larger than the cor- 
responding errors, we check our result with the 
X 2 -method Q described in the following. 
Let us review the general form of the scaling 
relations for different observables O 

O = LP/ u Q {tL^ v ,L- u ), (16) 

where we only take the largest irrelevant expo- 
nent into account. Here O is M , \ or g r with 
p = — /?, 7 and respectively. Expanding the 
function Qo to first order in the variables we 
hnd 

O = L^(co + (ci +c 2 L-")tL l ' v 

+ c 3 L-»), (17) 

which reduces to 

O = L p ' v (c + caL—) (18) 

at the critical point t = 0. Therefore at the 
critical coupling a fit with equation l|18[) has 
the minimal x 2 - Since we do not know the in- 
fluence of csL~~ u ' we started without this cor- 
rection term leaving L = 12 out. Fig. |3| (a) 
shows the result, we find a deviation from our 
preliminary result in case of the magnetisation 
M and the susceptibility Xv ■ The minima from 
the Binder cumulant and the susceptibility \ 
however coincide at J w 1.42865. 
We thereafter made fits with the correction 
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FIG. 3: The minimal per degree of freedom 
for fits according to equation 1181 1 for M (red), 
Xv (pink), x (blue) and g T (green). Figure (a) 
shows the results without a correcting term oc , 
whereas figure (b) uses u> = 0.5. The dotted black 
lines show our previous value of J c , the arrows the 
result of [l^| with its error bars. The curves of M 
and Xv lie on top of each other. 



term in the range lo = 0.5 — 1.5. The mini- 
mal x 2 an d a perfect agreement of J c for all 
observables is found at uj = 0.5. This result is 
plotted in Fig. 0(b). x 2 /d.o.f. increases with 
lo and shifts J c in case of M and Xv to smaller 
couplings, while the position calculated from x 
increases and J c from g r remains nearly con- 
stant. Since the fits get worse we can exclude 
w-values larger than 0.8. The positions of J c 
for 0.5 < lo < 0.8 coincide within the error bars 
ofEq. (T5). 

At the critical point the Binder cumulant has 
the form 

g r (L) = g r (J c ) + c 3 L-», (19) 
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FIG. 4: The magnetisation M as a function of the 
lattice extension L at the critical point T c and H = 
0. The dashed line in is a fit to the ansatz l|21|l with 
lu = 0.5. 
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FIG. 5: The derivative of the binder cumulant 
dg r /dJ as a function of the lattice extension L at 
the critical point T c and H = 0. The dashed line 
is a fit to the ansatz 124H with dz — since the 
corrections are negligible. 

with the universal value g r {Jc) and a small cor- 
rection term c 3 L~ u . For fits with different u> 
we find g r 

g r (J c ) = -1.94456(10). (20) 

The quality of the fits does not change much 
(x 2 /d.o.f. w 0.2 - 0.3) with different w so a 
better estimate of uu is still not possible. 

B. The critical exponents 

Since we now know the critical coupling we can 
study the finite-size behaviour of several ob- 



servables with Eq. I|I8|) . These can be extracted 
from our reweighted data at T c . The studied 
scaling relations are 

M = L-V"(ao + a 3 L-*), (21) 
X = V' v (bo + b 3 L-") (22) 

with u> € [0.5; 1.0[ and the exponents —fi/v and 
j/u as free parameters. Since these two ratios 
are connected by the hyperscaling relation 

1 = 3 - 2 A (23) 
v v 

it is necessary to study a further observable, for 
example the derivative of g r , which is given at 
T c by 

^ = L^(d + d 3 L-n- (24) 

fn this way two independent exponents (e.g. (3 
and v) can be estimated. We fit all observables 
in the range L = 12-72. From Eq. (|2"T|) we 
obtain 

- = 0.519(2), (25) 
v 

in which the error also includes an w-variation 
in [0.5, 1.0[. Here a larger ui shifts f3/v to a 
smaller value at nearly constant x 2 /d.o.f. ~ 
0.4. Fig. 0] shows the result with uj — 0.5. 
Our x _ fits yield 

- = 1.961(3). (26) 

Our results of (3jv and 7/z/ are tested with the 
hyperscaling relation 

2^ + 1 = d, (27) 
v v 

with d = 3 being the dimension of the model. 
The left hand side of this equation is 2.999(5), 
correct within the error. 

Finally we analyse the derivative dg r /dJ of the 
Binder cumulant at T c . This observable is di- 
rectly calculated from the spline connection of 
our reweighted data in the neighbourhood of 
the critical point. The errors are obtained with 
the jackknife method, which seems to underes- 
timate the errors, so we therefore assume the 
largest error of the different L-values for each 
lattice. Our fit to ansatz (|24l) without correc- 
tions to scaling (x 2 / d.o.f. k 0.3) is shown in 
Fig. For d 3 — we find 

- = 1.223(5). (28) 
v 



G 



Source 






exponent 


this work 


[12] 


[15] 


[16,-17] 


dg r /dJ 


l/v 


1.223(5) 


V 


0.818(5) 


0.819(3) 


0.790 


0.819 


M 


p/v 


0.519(2) 


P 


0.425(2) 


0.424(5) 


0.407 


0.424 


X 


l/v 


1.961(3) 


7 


1.604(6) 


1.608(4) 


1.556 


1.609 



TABLE II: The critical exponents for the 0(6) model estimated in this work compared to the theoretical 
work of Butera and Comi |l2|. Antonenko and Sokolov and Gracey [H[l7|. 



The final results of the critical exponents are 
summarized in Table ITT1 ft and 7 are calculated 
with the result of v and the ratios l|25[) and (|26() . 
The three last columns of the table show the 
results from ^ and [lSE^. Butera and 
Comi as well as Gracey are in good agreement 
with our values, but the results of Antonenko 
and Sokolov are farther away. 
In the following Sections we use the fixed crit- 
ical exponents /3 = 0.425 and v = 0.818. The 
remaining critical exponents arc calculated by 
the respective hyperscaling relations between 
the critical exponents. For u we will use the 
value u> = 0.5, which seems to be the best esti- 
mate in all investigations. 



III. SIMULATIONS AT H > 

The magnetisation M is now calculated from 
equation A transversal and a longitudinal 
susceptibility can be defined as 

V({ M 2 ) - M 2 ) 

= V( {^ f ) . 



XL 
XT 



(29) 
(30) 

We simulated at several constant J-values and 
increasing magnetic field, starting at H = 
0.00025. The used lattice sizes were L = 
24, 36, 48, 72, 96 and 120. Around 20, 000 mea- 
surements were performed in the (J, i?)-regions 
we used for our fits. The only exception was the 
data of L = 120, where we performed 10, 000 
measurements at J c and 5, 000 measurements 
at all other J-values. The integrated autocor- 
relation time for the energy and the magnetisa- 
tion is strongly dependent on the used J-values. 
At J c and J > J c we increased the number of 
cluster updates between two measurements to 
have autocorrelation times Ti„ t <6. 
In the symmetric phase (J < J c ) the situ- 
ation is different. While the measurements 
of the magnetisation are less correlated with 
Ti„t(M)<4, the correlation of the energy in- 
creases rapidly with decreasing H and J. It 
reaches values of Ti n t(E) <30 for the larger lat- 
tices. 



A. The critical isotherm 

At the critical point the critical scaling of the 
magnetisation is given by 

M(T e ,H) = d c H^ 5 (l + d\H uu °) , (31) 

where non-analytic corrections from the leading 
irrelevant scaling field are taken into account. 
They are not negligible in our model. The crit- 
ical exponents 8 and v c are known from the 
hyperscaling relations and only depend on the 
ratio P/v = 0.519(2): 

5 = 3~ - 1 = 4.780(22), (32) 
P 
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0.4031(24). 



(33) 



In order to exclude finite-size effects we carry 
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FIG. 6: The magnetisation at the critical coupling 
as a function of H. The solid line is the fit to the 
ansatz 131II . while the dashed line is the leading 
term. 

out a reweighting analysis for all lattices and 
fit the result from the largest lattice to approx- 
imate the value of V — > 00 . This is done for the 
interval H G [0.00075; 0.04] and we find 



0.642(1). 



(34) 



7 



2.0 - 



1.5- 



1 . 



. 5 



. 



I I I I 1 


1 


1 




" MT ^ ' y 









J? 

X 

* 




: 12 

+ 96 ~ 
72 

x 48 ; 





" « 

' & 




36 

o 24 













- s 








" 









o 






hl 17 ^: 





50 100 150 200 

FIG. 7: The finite-size-scaling function Qm,x ° n 
the critical line, Eq. 1371 , The solid red line indi- 
cates the asymptotic function Qo.oo for z'>40. 



Our result is plotted in Fig. El There are min- 
imal negative corrections. If one treats 5 as 
a free parameter the result 6 — 4.79(1) agrees 
with our first estimate. 

The finite-size-scaling function for the magneti- 
sation is 

M(T,H,L) = L~ fJ/u ■ 

<5>{tL 1/ \HL 1/y %L- bJ ).{?,h) 

The scaling function $ can be expanded in L~ u 
to 

M(T,H,L) = L-^ v ® {tL l l v ,HL l ' u ') 

+ ... . (36) 

At T c the leading part is now given by 

M(T C ,H,L) = L-^Qm(z') (37) 

with the universal scaling function Qm{z') and 
the argument z' = HL 1 ^". The results of all 
lattices are shown in Fig. The data points 
scale very well and the influence of corrections 
to scaling is small. In the limit z' — > oo we 
expect the asymptotic behaviour 

QmMz') = d c z n/s , (38) 

which is observable for z'>40. This way one 
checks the value of d c with a fit of reweighted 
z'-data of the larger lattices L — 72, 96, 120. 
We find d c — 0.642(1), which agrees perfectly 
with our first value of cL. 



B. Numerical results at T / T c 

Let us review some perturbative predictions for 
the magnetisation and the susceptibilities. The 
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FIG. 8: The magnetisation in the broken phase 
as a function of \[H for the couplings J = 1.6, 
1.55, 1.50, 1.47, 1.45 and J c and different lattices, 
starting with the largest J-value at the top. The 
solid lines represent interpolations from a reweight- 
ing analysis of the data. The dashed lines are the 
fits to ansatz 141H while the dotted line is the fit of 
equation (I31H at J c . 



continuous 0(6) symmetry of our spin model 
gives rise to spin waves, which are slowly vary- 
ing (long-wavelength) spin configurations with 
energies arbitrarily close to the ground-state 
energy. In d > 2 they are massless Gold- 
stone modes associated with the spontaneous 
breaking of the rotational symmetry for tem- 
peratures below the critical value T c [18|, For 
T < T c the system is in a broken phase, i.e. the 
magnetisation M(T, H) attains a finite value 
M(T, 0) at H = 0. 

The transverse susceptibility has the form 



XT 



M(T, H) 
H 



(39) 



for all H and T. This relation is a direct conse- 
quence of the 0(6) invariance of the zero-field 
free energy and can be derived as a Ward iden- 
tity [3. 

The longitudinal susceptibility diver ges on the 
coexistence curve for 2 < d < 4 pol 12 1|. The 
leading terms in the perturbative expansion for 
three dimensions are 



X l(T<T c ,H) = b (T)H 



-1/2. 



■c 2 (T). (40) 



Since the susceptibility is the derivative of the 
magnetisation with respect to H we find for the 
magnetisation 

M[T < T C .H) = M(T, 0) + Cl (T)H 1/2 

+ c 2 [T)H (41) 
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J=l/T 


M(T,0) ci(T) 


c 2 (T) 


10 4 ■ H 


X 2 /dof 


1.45 


0.1701(03) 1.339(15) 


-2.86(18) 


9-25 


0.45 


1.47 


0.2219(02) 0.924(02) 


-1.138(14) 


12-74 


0.56 


1.50 


0.2761(01) 0.659(01) 


-0.436(08) 


10-93 


0.27 


1.55 


0.3401(01) 0.463(01) 


-0.141(02) 


10-163 


0.78 


1.60 


0.3878(01) 0.363(01) 


-0.047(01) 


19-175 


0.49 



TABLE III: Parameters of the fit of M(T < T c ,0) to the ansatz igTJl. The fifth column is the used fit 
range. 
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FIG. 9: Magnetisation M(T,0) on the coexistence 
curve as a function of (T c — T). The red line is the 
fit to ansatz 14211 while the black dashed line is the 
leading part. 
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FIG. 10: The magnetisation in the symmetric 
phase (T > T c ) as a function of H, starting from 
the top with fixed J = J c , 1.42, 1.41, 1.40 and 1.38 
and different L- values. The lines are spline connec- 
tions between the data points. 



near the coexistence curve. Fig. 03 shows 
our results of the magnetisation in the bro- 
ken phase and the corresponding extrapola- 
tions to M(T, 0) in the thermodynamic limit 
(V — > oo). The numbers of the parameters 
are presented in Table ITTT1 The ii-extension of 
the regions, where the predicted Goldstone be- 
haviour is found, increases with J, while finite- 
size effects become larger at small H and larger 
J (L>160 would be necessary for finite-size in- 
dependent data). 

We fitted the values of M (T, 0) to the form 

M(T < T c , 0) = B (T c - Tf[l + b x (T c - Ty v 
+ h (T c - T)] (42) 

with fixed values j3 = 0.425, ujv = 0.409 and 
the result 

B = 1.22(1), (43) 

bi = -0.184(49) and b 2 = 0.31(13). The error 
of B also includes the slight uncertainty in the 
value of uiv. Our final result of M(T, 0) and 



the difference to the leading term are plotted 
in Fig. El 

Since one of the main aims of this work is the 
determination of the magnetic equation of state 
in Section IIVI we also simulated in the high 
temperature phase. Again we use data of the 
largest lattices as an approximation for the infi- 
nite volume value. The result is plotted in Fig. 
1101 From a Taylor expansion we expect 

M(T >T C ,H) oc H. (44) 

at small H . With increasing temperature the 
H- interval with this behaviour increases, as one 
can see in Fig. ^| 

IV. THE SCALING FUNCTION 

The critical behaviour of the magnetisation in 
the vicinity of T c is described by the general 
Widom-Griffiths form |23 

y = f(x) (45) 
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FIG. 11: The function M/h 1/s . The solid lines 
in the broken phase are the reweighted results for 
M/h 1/s at J = 1.60, 1.55, 1.50, 1.47, 1.45, 1.445 and 
1.44, from the bottom to the top. They are ex- 
trapolated with Equation I55|l to fa (dashed line). 
The circles are single data points in the symmetric 
phase. 



with 



y = h/M s , x = t'/M 



V/3 



(46) 



where the variable t' is proportional to (T — T c ) 
and h proportional to H. A common normali- 
sation of the function f(x) is 



/(0) = 1, /(-i)=o. 



(47) 



The variables t' and h are the conveniently nor- 
malized reduced temperature t' — (T — T c )/Tq 
with T = B- 1 /' 3 = 0.626(12) and the reduced 
magnetic field h = H/Hq using Hq = d~ s = 
8.32(6). The function f(x) is universal and was 
derived from the e-expansion (e = 4 — d) to or- 
der e 2 19]. In the limit x ->• -1, i.e. at T < T c 
and close to H = the result was inverted to 
give x + 1 as a double expansion in powers of y 
and y d ' 2 - x H2 



1 = c lV + Z 2 y d l 2 - X + d iy 2 
+ d2V d/2 + d 3 y d - 2 + 



(48) 



The coefficients ci , C2 and are thereafter ob- 
tained from the general expression of |19j . 
In the large-x limit (corresponding to T > T c 
and small H), the expected behaviour is given 
by Griffiths's analyticity condition 



fix) 



„ 7 -2(n-l)/3 



(49) 



The form (|45|l of the equation of state is equiv- 
alent to the often used relation 



M = hV'fdz), 



(50) 



where fa is a further universal scaling function 
and z the combination 

2 = t'/h 1 ^ 5 . (51) 

The normalisation conditions of fo( z ) are 



/ G (0) = 1 and f G (z) 



(-zf . (52) 



This version is normally used for comparison 
to QCD lattice data. The function f(x) is con- 
nected with fdz) by 



y = f G s > 



x = zf 



-1//3 



G 



(53) 



These scaling functions are only valid close to 
J c and H — > 0. First tests show that the data 
we have used in the broken phase does not scale 
directly, while in the high temperature phase 
most of the data scales close to T c and small 
H . So we used a more general form of l|50|l 



Mh~ 1/S = ^(z,h uv ") 



(54) 



with a scaling function ^P, which can be ex- 
panded to 



Mh- 1/S = f G (z) + hP^fgiz) 



(55) 



This way in the broken phase one obtains the 
leading part fa by quadratic fits to our data in 
h UVc at constant z-values and different (J/H)- 
combinations. But we are only able to correct 
the data with z < —2 because we have not 
enough J-values closer to J c to make the fits. 
In Fig. ^] we show the influence of the correc- 
tions and the final scaling function Jq in the 
broken phase (dashed line). 
Our result for fciz) can be transformed with 
(I53|) into the Widom- Griffiths form of the equa- 
tion of the state H45|) . Unfortunately the z- 
interval we used to extract fc(z) is only equiv- 
alent to the small region —0.9 < x < —0.7, 
which can be used for a fit. We use the three 
leading terms in l|48|) 

xi(y) + 1 - (ci + d 3 )y + c 2 y^ 2 

+ d 2 y 3 / 2 . (56) 

Since y(0) = 1 the coefficients are connected by 
d 2 = 1 — (ci + d 3 + c 2 ). Fits to x in the interval 
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FIG. 12: The functions y = f(x) at small x (a) and 
large a;- values (b). In figure (a) we plot the fit to 
ansatz l|56jl using extrapolated data in the broken 
phase and data-points in the symmetric phase. In 
(b) we plot data-points of the symmetric phase and 
the fit to ansatz 14911 using the first three terms. 
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FIG. 13: (a) The scaling function f G of the 0(6) 
model (green) and the asymptotic behaviours for 
z — > oo (black) and z — * — oo (red), (b) The scaling 
function fa for the Ising (black), 0(2) (red), 0(4) 
(blue) and 0(6) model (green). 



—0.9 < x < —0.7 and points in the symmetric 
phase with 0.2 < x < 2.9, 1.4 < J < J c and 
H < 0.0015 lead to 

c x +d 3 = 0.36(5) , c 2 = 0.69(3). (57) 

The result of the fit is shown by the line in Fig. 

For large x we use a 3-parameter fit of the first 
three terms of Griffiths 's analyticity condition 
(£21 

y 2 (x) = a lX ~< + a 2 x^- 2p + a 3 x^ il3 (58) 

in the interval x £ [1.75,202] and data points 
restricted to 1.4 < J < J c and H < 0.0015. We 
find 

ai = 0.92(1), a 2 = 1.17(2), 

a 3 = 0.91(3). (59) 



This result is plotted in Fig. U^b). With the 
coefficient ai of the leading part one can calcu- 
late the universal ratio 

R x = a^ 1 = 1.09(1). (60) 

The 0(6) scaling function fa can be paramet- 
rically obtained from x\(y) and y 2 {x), which 
is connected by a spline between z — 0.5 and 
z = 0.8, where we have no reliable parametri- 
sation. The result is plotted in Fig. ll3f a). Also 
plotted are the leading terms of the asymptotic 
behaviour at z — > ±oo. These are 

fa(z) z ^°° {-zf , (61) 
f G (z) 2 ^°° R x z-^ (62) 

according to the normalisation (|52|l . The fact 
that for large temperatures and small H the 
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FIG. 14: The scaling function f x for the 0(6) 
model. The dashed red line shows the position of 
z n = 1.34. 



magnetisation is proportional to H, see Eq. 
(|44|l . explains the asymptotic behaviour for 
z — > oo. In the symmetric phase the asymptotic 
behaviour is reached for small absolute values 
of z, while in the broken phase the scaling func- 
tion converges to the asymptotic form not until 
large absolute values of z. 

Finally, the 0(6) scaling function fa can be 
compared to the corresponding functions for 
the Ising (0(1)) model H|, the 0(2) [IcJ and 
the 0(4) model 0, shown in Fig.EJb). All 
functions have a similar shape. 



THE PSEUDOCRITICAL LINE 
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pseudocritical line, (b) is a double-log plot of (a). 
The solid line in (b) shows the asymptotic form 
Qz p , oo, the symbols denote different lattice sizes L. 



In order to discuss finite-size-scaling functions 
in an easier way it is common to study lines of 
constant z- values. There one expresses H as a 
function of T or vice versa. Important examples 
of lines of fixed z are the critical line [z = 0), 
discussed in Section ITlI Al and the pseudocrit- 
ical line z = z p = const, the line of peak posi- 
tions of the susceptibility xl in the (i, /i)-planc 
for V — y oo. There are two different ways to 
find that value of z p for O(N) models. One 
way is to locate the peak positions of xl as a 
function of the temperature at different fixed 
small values of the magnetic field on lattices 
with increasing size I? . The scaling function, 
on the other hand, offers a more elegant way to 
determine the pseudocritical line. Since xl is 
the derivative of M 



its scaling function f x (z) can be calculated di- 
rectly from fc{z) 

/xW = i(/oW-|^)) ■ (64) 

The maximum of f x (z) is located at z p , which 
is another universal quantity. We find for the 
0(6) model 

z p = 1.34(5). (65) 

The error includes the fit-errors of the parame- 
ters in 1 (55) 1 . 

In Fig.^Jwe show the result for f x (z) from Eq. 

in the 0(6) model. 
At finite-size-scaling analysis in the vari- 
ables H and L can be performed. Eq. (|36J) 
reduces to 

M(H,L) = L-V" Q z (hL x l v °) + ... (66) 
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with another universal scaling function Q Zp . 
The asymptotic form Q 2p ,oc of Q Zp is 

Q Zp L =^ Q Zp ,oo = !g{z p ) (hL 1 ^) 1 / 5 (67) 

The results are presented in Fig. I15f a). The 
data does not scale directly but with increas- 
ing volume the data-points approach Q Zp from 
the top. In Fig. HaT h) we plotted the data 
in a double-log form and find that the asymp- 
totic form Q Zpt oo coincides with the Q 2p -value 
of the largest lattice extension at HL 1 ^" > 42. 
Therefore Q z is asymptotic. At smaller values, 
one observes an approach of Q Zp from below to 

Qz p , oo- 

VI. CONCLUSIONS 

In this paper we calculated several important 
quantities of the 0(6) spin model directly from 
Monte Carlo simulations on cubic lattices. At 
zero external field we determined the critical 
coupling J c by a finite-size-scaling analysis of 
the Binder cumulant and by the x 2_me thod. 
Our result agrees in the first four digits with 
the result of Butera and Comi. At the critical 
point we estimated the critical exponents from 
finite-size-scaling fits. We obtained f3 from the 
magnetisation, 7 from the susceptibility and 
v from the derivative of the Binder cumulant. 
Our results are in accord with the values found 
by Butera and Comi but slightly different 
compared to the values found by Antoncnko 
and Sokolov. We find small corrections to 
scaling for all observables. 
On the critical line T = T c , H > and in the 
limit V — > 00, the critical amplitude d c of the 
magnetisation was computed. We found small 
negative corrections to scaling and checked the 
finite-size-scaling behaviour of M at T c and its 
asymptotic form. 

Below the critical temperature, we investigated 
the behaviour of M at several couplings J as 
a function of H 1 / 2 in the limit V — > 00. Close 
to the coexistence line, i.e. small H — * 0, the 
predicted Goldstone behaviour was observed. 
We were able to extrapolate our data to the 



values M(T < T c ; H = 0) of the infinite volume 
limit, fitted these M-values with the corre- 
sponding ansatz, and estimated the critical 
amplitude B of the magnetisation. In this case 
the corrections to scaling were again negative 
and more pronounced as on the critical line. 
At high temperatures and H > 0, we observe 
the expected proportional dependence on H of 
the magnetisation. 

We used our data of the largest lattices in the 
low and high-temperature phase to parametrise 
the scaling function /g of the 0(6) model. We 
encountered large corrections to scaling in the 
broken phase, while most data in the symmet- 
ric phase scales directly. By generalizing fc 
to include corrections to scaling, our group 
extracted a part of fa in the broken phase 
and fitted the result combined with direct 
data points in the symmetric phase. On the 
other hand, we fitted data of the symmetric 
phase using Griffiths 's analyticity condition. 
Finally, we compared our 0(6)-result for fc 
with the corresponding scaling functions of the 
0(1), 0(2) and 0(4) model. These functions 
are clearly distinguishable and in a systematic 
order. We use our result of fc to calculate the 
scaling function f x of the susceptibility. From 
the position of the maximum in f x the location 
of the pseudocritical line was determined. 
There we made finite-size-scaling plots and 
found considerable corrections to scaling. The 
data of smaller lattices approaches the uni- 
versal finite-size-scaling function from above. 
The asymptotic form of the universal part is 
reached at HL l / Vc « 42. 

A comparison between the 0(6) model and 
aQCD will be done in the near future. 
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